Transverse Oscillating Waves and the Effects of Capacitive Coupling in Long Intrinsic 

Josephson Junctions 



Farzad Mahfouzi 

Department of Physics and Astronomy, University of Delaware, Newark, DE 19716-2570, USA 

In this paper we investigate the excitation of longitudinal and transverse plasma waves in intrinsic 
Josephson junctions. We consider the outermost branch of IV characteristic (IVc) in current biased 
case and try to find the conditions in which plasma waves can be excited. We change the parameters 
of the system and get the corresponding breakpoint current at which the plasma waves start to 
initiate. We present specifically the modes containing only transverse waves where we can have 
radiation. As a result we find the range of parameters that the system can radiate. 
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One of the topics in high Tc superconductors that 
has recently attracted a great interests in the litera- 
ture is the observed terahertz wave emission from the 
layered structured high-temperature superconductors, [ij 
It's widely believed that cuprate high-temperature su- 
perconductor hke Bi2Sr2CaCu208 (BSCCO) has a lay- 
ered structure which leads to an intrinsic Josephson junc- 
tion (IJJ) effect. Applying current in z-direction we 
get AC Josephson oscillations that works in terahertz 
range of of frequencies. Moreover, the parameters in in- 
trinsic Josephson junctions are controlled by the atomic 
crystal structure rather than by leads and amorphous 
dielectric layer in artificial JJs. This results in the ex- 
istence of different mechanisms of coupling among IJJs 
e.g. phonon coupling, capacitive coupling,^ inductive 
coupling, charge imbalance effect [4j]. In this paper we 
consider inductive-capacitive coupling model where we 
can have longitudinal and transverse oscillations. Radi- 
ation from this system has already been investigated in 
literature d, 0| ■ In this paper we investigate properties 
of the transverse and longitudinal plasma waves excited 
by a constant bias current in the multi-Josephson junc- 
tions. We restrict ourselves to the two-dimensional IJJs, 
for simplicity, and consider that the length of the sys- 
tem in y-direction is small comparing to the magnetic 
penetration depth parallel to the ab-plane of the crystal. 
Therefore we assume the system to be homogeneous in 
y-direction. The layers are stacked along the z-direction 
(c-axis). We consider that the length of the junctions in 
x-direction is comparable to the magnetic penetration 
depth and the phase-differences of junctions along the x- 
direction varies accordingly. It has been established that 
the junctions in these multi-junction systems are coupled 
both capacitively and inductively with each other. We 
show that the plasma waves inside the junctions are ex- 
cited parametrically by the Josephson oscillations arising 
from the bias current. 

The couplings can be incorporated into the dynam- 
ics of the phase-differences in IJJs, using the generalized 
Josephson relations. 0,0,11^ = x (2^) " 
acV<^^^Vy,^u.,){x)) and ^ = i^(Hy[^^[^_^^{x) - 



ah'^'^'^^ Hy ^i i_i){x)). Where the parameters ac and ul 

are respectively the capacitive and inductive coupling 
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constants given by ac = ^ , Q^l = -jj- with s,d,£,fi 
andAab being, respectively, the thickness of the supercon- 
ducting layer and insulating layer, the dielectric constant 
of the insulating layers, the charge screening length of 
the superconducting layers and the magnetic penetration 
depth parallel to the ab-plane of the crystal. The equa- 
tion of motion for the gauge-invariant phase-differences 
ip{x,t) in dimension-less form is 0,0,0] 

X^iL~%^ - iC~%^ = sin(^.) + /3v3. - I/L 

(1) 

where C = l-acV^^) and L = I-olV^^), = ^ejlyi, ' 
T — u!„t and uj?, — ^^Hr- The index indicates the ith 
junction and V'-^^' is the abbreviation of the second rank 
difference. This equation has been obtained from con- 
servation relation where the total dimension-less current, 
I / Ic, can be written in terms of the transverse current, 
T displacement current, ^ — Ohmic cur- 
rent, (3ipi, and the Josephson current, sin^ipi). 

In this model we consider that an electron in the posi- 
tion X of ith layer can only jump to the same position on 
the next layer i.e. J^^i+i — sin{ipi) + f3ipi. We call this as 
local hopping approximation which is an appropriate ap- 
proximation in cases of slow variation of phase difference 
in X direction. In this paper we use this approximation 
and neglect the fast varying solutions. 

In the following we are going to get some information 
about the system of equations (1). It's well known that 
in the case of short Josephson junctions the IVc has a 
multi branch structure corresponding to different config- 
urations of the rotating and oscillating states among the 
junctions in z direction. It has been shown that neglect- 
ing the effects of boundary layers, the number of branches 
is equal to the number of the Josephson junctions. [1| But 
what about the long Josephson junction case? Does it 
lead to the same number of branches? In order to answer 
this question, let us first solve the equation of motion for 
one Josephson junction. The result of this solution is 
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FIG. 1: (Color online) IV characteristic of one Josephson 
junction which is equivalent to the case q:l = 0, ac = 0. 
As the left inset of the figure, we show the frequency spec- 
trum of the derivative of voltage with respect to x and as the 
right inset of the figure we show the distribution of voltage 
along the junction at an arbitrary time in the regions. The 
Current is defined in units of Josephson critical current and 
the Voltage is defined in units of huip/e. 



presented in Fig. ([T|). To obtain this result we have 
considered 100 points along the junction {x direction). 
Then using fourth order Rounge-Kuta algorithm we have 
solved the system of equations in time. This result shows 
another kind of branch structure which is well known and 
corresponds to different distributions of phase along the 
junction. In general we can distinguish two types of solu- 
tions corresponding to flux-flow mode and oscillating dis- 
tributions. It's clear from this result that if we start from 
outermost branch we reach first the oscillating mode and 
then to the flux-flow mode. To clarify the situation we 
specify some regions of the IVc by letters (A), {B), (C), 
(D), (E), (F) and present, as left inset, the FFT analy- 
sis of the time dependence of ^('^-O-O^I^'^C^-Oi) units of 
Josephson frequency and as right inset, the distribution 
of voltage along the junction. We can see three, two and 
one flux-flow mode states at the regions (C), (D) and 
{F) respectively and standing oscillating waves at the re- 
gions {B) and {D). Although the number of fluxons can 
be seen directly from the distribution of voltage, we can 
also find it from the number of peaks in unit range of the 



FFT results i.e. N^,^ 



Let us now consider the multi layered structure and 
find the number of branches in IVc in cases of Long 
Josephson junctions. From above discussion we find that 
there are at least two mechanisms of branching in IVc 
j.e.(l) different configurations of the rotating and oscil- 
lating states among the junctions in z direction and (2) 
different distributions of phase along the junction in x 
direction. Should we consider both effects the number 
branches exceeds the number of junctions which is ap- 



parently not in accordance with experimental results. In 
the following we see that in the case of high numbers of 
junctions the number of branches tends to the number 
of junctions. We show this by noting that at large num- 
ber of junctions, upon creation of longitudinal plasma 
wave at the breakpoint current the system immediately 
transits to another branch containing less resistive junc- 
tions. Therefore, should the system starts from an os- 
cillating mode with longitudinal waves, the system does 
not reach other resonance regions corresponding to some 
other transverse modes and so the additional branches 
due to inductive coupling in x-direction does not appear. 
But what if the system in the outermost branch starts 
from a state with only transverse mode (without longi- 
tudinal mode)i.e. (n,0)? It's clear that, in this case, upon 
decreasing the current, the system continues to be in this 
state until it reaches to another resonance region which 
contains longitudinal plasma wave and then the system 
jumps to another branch corresponding to less number 
of resistive states. 

From Fig. ^ we learned that the flux-flow modes ap- 
pear far from the breakpoint current. Additionally we 
know that in multi Josephson junction case, upon cre- 
ation of longitudinal plasma waves the system jumps to 
another state with less number of resistive junctions. In 
the following, we neglect the flux-flow modes and con- 
sider only the region close to breakpoint current of out- 
ermost branch in IVc. Later in this paper we can see 
that, if we consider the longitudinal coupling, the system 
in fact hardly reaches to the flux-flow regions and there- 
fore the flux-flow modes play no effective role in multi 
Josephson junctions. As a result we start from points 
(A) and (B) in Fig. ([T]) and investigate their character- 
istics. Point (A) represents breakpoint current and the 
system at this current enters into the parametric reso- 
nance region. Throughout this paper we consider this 
point only and discuss about its features. In Fig. ^ we 
see that, the region (B) which we call radiation branch, 
has close characteristics to the point (A) and also it ap- 
pears after the breakpoint current. In the following we 
discuss that this branch might be observable in IVc. 

By the fact that on the outermost branch all of the 
junctions are in the resistive (rotating) state, we neglect 
the effect of boundaries (large number of junctions) and 
consider that initially, phases of the junctions are the 
same. We consider also that this phase does not depend 
on the position x i.e. ipi{x,t) = 4'{t). Then we test the 
stability of such a state by adding a small term to the 
phase i.e. ipi{x,t) — (j)(t) +di{x,t). Linearizing the equa- 
tion with respect to this term and then making Fourier 
transformation with respect to x and i leads 
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FIG. 2: (Color online)/? dependence of the breakpoint current 
in units of critical current ai ul ~ l,e — 1 and different 
capacitive coupling (ac) 



|l + ^| + ^i«W = W (3) 

where C(fc^) ^ l + 2ac(l-cos(fc^)), L{k^) = l + 2aL(l- 
003(^2)), e — ttA^/L^ and n is an integer number. Due 
to absence of external magnetic field the derivative of 
phase at the boundaries of the junctions are zero. This 
boundary condition leads to k,j; = mr. On the other hand, 
considering a periodic boundary condition for z-direction 
we get kz = 2m'K /N . Where m is also an integer number 
and N is the number of junctions, n and m represent the 
transverse and longitudinal wave numbers respectively. 

We find that the this initial condition is stable at the 
big currents and unstable bellow some critical current 
which we call breakpoint current. We find the breakpoint 
current at different parameters of the system. Fig. [2] 
shows breakpoint current versus (3 at different capacitive 
coupling parameters. The black curve shows the max- 
imum of breakpoint current for different modes (n, m) 
where n and m arc taking the values 0,1,2. Where 
fnmax = 2 means that the number of Josephson junctions 
is 2 and rimax — 2 means that we are neglecting fast vary- 
ing solutions. Additionally one can see that bigger values 
of n appear at smaller /3. Therefore, to see the results 
corresponding to bigger n we need to choose smaller /3 
which is out of scope of this paper, because, we are ne- 
glecting fast varying solutions. In Figs. [21 [3l [4] we show 
especially the modes that longitudinal waves are absent 
i.e. the oscillation that can lead to radiation. We note 
that the parts of these curves that coincide with the black 
curve are observable only. We note also that increasing 
the number of junctions smears the /3-dependence of the 
breakpoint current. 

From experimental point of view, it might be interest- 




FIG. 3: (Color online)/? dependence of the breakpoint current 
in units of critical current at ac = l,e = 1 and different 
inductive coupling (ai) 




FIG. 4: (Color online)/? dependence of the breakpoint current 
in units of critical current at etc = l,aL = 1 and different e 



ing to know the width of the region in IVc that there's 
electromagnetic radiation and it's dependence on the pa- 
rameters of the system. The resonance width is defined 
by the difference between breakpoint current and jump 
point current at which the system goes to another branch. 
We have already shown that in the presence of longitu- 
dinal plasma waves the width of the resonance region is 
getting zero with increasing the number of junctions, [l^ 
In case of transverse plasma waves the situation is differ- 
ent. In addition to the fact that the breakpoint current 
in this case doesn't depend on the number of junctions, 
result of simulation for one Josephson junction shows also 
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FIG. 5: (Color online)Width of radiation range on the outer- 
most branch of IVc {l^I Radiation) versus P at different capac- 
itive coupling etc 



FIG. 6; (Color online) Minimum frequency of the plasma wave 
spectrum in units of Josephson plasma frequency ( black 
curves) and breakpoint Josephson frequency (red curves). 



that the return current doesn't depend strongly on the 
mode of transverse plasma waves. This leads to the con- 
clusion that jump point current can be obtained from the 
breakpoint of the modes with the longitudinal plasma 
waves. Therefore we can get the range of current that 
we have only transverse modes and and can be touched 
from the outermost branch. We can see the result of this 
calculation in Fig. [5l This model predicts that if we 
get the branch corresponding to radiation (only trans- 
verse waves) in the return current part of the outermost 
branch, we can then increase the current and have radi- 
ation until the end of the branch. But one should note 
that the end of branch in this case is not critical current 
[I — Ic). In addition, if we include the dissipation due 
to the radiation power, it will lead to increase in volt- 
age in IVc and therefore decrease in the length of the 
radiation branch. We note that although increasing the 
number of junctions has not strong effect on the peak 
of the /3-dependence of l^I Radiation but it decreases the 
range of (3 that Aluadiation is nonzero. This shows that 
increasing the number of junctions makes the situation 
more difficult to get the radiation region. 

Now that we know the condition for electromagnetic 
radiation one question arises; What is the frequency of 
oscillation at the breakpoint current and breakpoint re- 
gion (wjbp)? Results of our numerical calculations e.g. 
(A) point in Fig. [H shows that the Josephson frequency 
at the breakpoint in IVc follows the {M+^)ujjbp relation, 
where M is integer number. We note that this result cor- 
responds to the breakpoint of outermost branch which is 
not stationary. The stationary solution which can be ob- 
tained from the original nonlinear equation and we call 
radiation branch, contains frequencies equal to ^Mujjbp. 



In Fig. [n]we show the first frequency of the plasma waves 
(M = 0) in units of Josephson plasma frequency ( black 
curves) and in units of breakpoint Josephson frequency 
(red curves). We can see that the frequency of plasma 
waves is close to half of the Josephson frequency at the 
breakpoint. The deviation is due to the fact that the 
system is current biased. 

We conclude that upon decreasing the bias from the 
fully resistive state, the system enters into resonance be- 
tween the Josephson frequency and the plasma waves. 
As a result we have stationary plasma waves that can 
contain longitudinal and transverse waves. Should the 
longitudinal wave exist, some of the junctions immedi- 
ately switch back from the resistive into superconduct- 
ing states. In case of transverse waves without longitu- 
dinal one, an additional branch, that we call the radi- 
ation branch, appears and the system starts to radiate 
electromagnetic power. Upon decreasing current on the 
radiation branch the system finally enters into the res- 
onance regions corresponding to other modes of plasma 
waves that yields a transition to another state with less 
resistive states. Therefore in the long intrinsic Josephson 
junction systems where both longitudinal and transverse 
waves exist, the branching in IVc is mainly due to longi- 
tudinal waves and the only observable effect of transverse 
waves is accompanied by radiation. 
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